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Moving, merging and annihilating Dirac points are studied theoretically in the tight-binding model 
on honeycomb lattice with up-to third-nearest-neighbor hoppings. We obtain a rich phase diagram 
of the topological phase transitions in the parameter space of direction-dependent hoppings. We 
obtain the conditions for the three Dirac points to merge and for the tricritical points. We find 
that only very small third-nearest-neighbor hoppings are enough for the existence of the merging 
of three-Dirac-points and the tricritical points, if the system is sufficiently anisotropic. The density 
of states is obtained to be D{e) oc |e| 1//3 when three Dirac points merge, and D(c) oc e^ 4 at the 
tricritical points. It is possible to realize these topological phase transitions in the ultracold atoms 
on the optical lattice, strained monolayer graphene or strained bilayer graphene. 



PACS numbers: 73.22.Pr, 73.43.Nq, 71.10.Pm, 73.20.At 



I. INTRODUCTION 

Massless Dirac fermions are realized in various fields 
in condensed matter physics and they attract much in- 
terest recently, e.g. graphene 1-3 , the organic conduc- 
tor a-(BEDT-TTF) 2 l3 4,5 and iron-based superconduc- 
tor, BaFe2As2 6 . These materials have Dirac points and 
the linear band structure near the Fermi energy. 

The moving, merging and annihilating Dirac points 
due to the breaking of the rotational symmetry were 
studied 7 and observed recently in the ultracold atoms 
on a optical lattice 8 . Many authors have theoretically 
studied the merging Dirac points in a-(BEDT-TTF) 2 l3 9 , 
strained or twisted bilayer graphene 10-14 , and the hon- 
eycomb lattice with third- nearest- neighbor hoppings 15 . 
Most of the studies have used the effective model, which 
describe the energy near the Dirac points. 

The energy band of graphene is obtained by first- 
principles band calculations 16-18 , and the band structure 
near the Fermi energy can be reproduced by a simple 
tight-binding model 19 . If hoppings are only between the 
nearest sites and they are independent of the direction, 
there are two Dirac points at K and K' in the Brillouin 
zone. Reich et al. 2a have shown that up to third-nearest- 
neighbor hoppings are necessary to describe the band 
structure in all Brillouin zone. They have obtained that 
the nearest-neighbor, the next-nearest-neighbor and the 
third-nearest-neighbor hoppings are 2.79-2.97 eV, 0.073- 
0.68eV and 0.3-0.33 eV, respectively 20 . 

A finite gap is important for the application to nano 
devices. There are two routes to open the finite gap at the 
Dirac points. One route is via the the violation of the in- 
version symmetry. The violation of the inversion symme- 
try can be caused by the difference of the on-site poten- 
tials on A and B sublattices or the sublattice-dependent 
second-neighbor hoppings 21,22 . The finite gap observed 
experimentally in graphene on a SiC substrate 23 is caused 



by the breaking of the inversion symmetry. 

The other route to open a finite gap is realized by 
making two Dirac points with opposite topological num- 
ber merge and annihilate 7,9-15 . When the hoppings be- 
tween nearest sites become different in three directions, 
the Dirac points move from K and K' to one of the three 
inequivalent M points (Mi, M2 and M3), merge into a 
scmi-Dirac point and annihilate to make a finite energy 
gap. In order to make a finite gap in monolayer graphene, 
we need the deformations of the order of 20% caused by 
the strong uniaxial strain 24 or the shear strain 25 , which 
are difficult to realize. However, the parameters can 
be controlled in the ultracold atoms on the optical lat- 
tice and the merging and annihilating Dirac points are 
observed 8 . 

It has been shown that the strained or twisted bi- 
layer graphene is modeled by the single-layer honey- 
comb lattice with the third-nearest-neighbor hoppings, 
the strength of which are near the half of the nearest- 
neighbor hoppings 14 . 

In this paper we study the tight-binding model on 
the honeycomb lattice with up-to third-nearest-neighbor 
hoppings, and we study the moving and the merging of 
the Dirac points by changing the strength of direction- 
dependent hoppings. We obtain interesting phase dia- 
grams, in which we can see the merging of three Dirac 
points and the tricritical points. The density of states 
due to Dirac points and merged Dirac points are calcu- 
lated. 

The model is given in Section II. The topological num- 
ber and the Berry phase of the Dirac points are discussed 
in Section III. In Section IV we study the phase diagram 
in which we consider only the phase transition due to 
the merging of two Dirac points at Mi, M 2 , M3, and 
I\ Moving Dirac points are studied in the model with 
direction-dependent nearest-neighbor hoppings in Sec- 
tion V. Rich phase diagrams in the model with direction- 
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FIG. 1. (Color online) Honeycomb lattice and hoppings. 
Filled circles (•) and open circles (o) are sites in A sublat- 
tice and B sublattice, respectively, vi and V2 are the unit 
vectors. The hoppings between the nearest sites, the next- 
nearest sites, and the third-nearest sites are shown in (a), 
(b), and (c), respectively, and they are considered to depend 
on the direction and the sublattices. 



dependent nearest-neighbor and third-nearest-neighbor 
hoppings are studied in Section VI. Section VII is the 
conclusions and detailed calculations are given in Ap- 
pendix. 



II. TIGHT-BINDING MODEL ON 
HONEYCOMB LATTICE 

We study the tight-binding electrons on the honey- 
comb lattice (see Fig. 1). There are two sublattices, A 
and B (open and filled circles in Fig. 1), which form tri- 
angular lattices. We take the unit vectors as 



i 



-a, 



(1) 



v 2 



73 f N 



and the vectors connecting the nearest sites as 
8 a = ~(v a + v 2 ) = (^a,0), 
<*& = g(v 2 - 2vi) = ( — —a, -a), 



and 



s c = ~(vi - 2v 2 ) = (-^o, -io), 



(2) 

(3) 
(4) 

(5) 



where a is the lattice constant. Hereafter, we set a = 1. 

We consider the direction-dependent hoppings between 
nearest sites (t a , t b , and t c ), the next-nearest sites (t' a , t' b , 
t' c , t", t'l and t") and the third- nearest sites (t 3a , t 3b , and 
t 3c ), as shown in Fig. 1. The hoppings between the next- 
nearest sites within the A sublattice are labeled by t' a , t' b 
and t' c , and those within the B sublattice are labeled by 
t^, t b and t" c . The Hamiltonian is given by 



where 



Ck = 



«k 
b k 



(6) 



(7) 



ak and 6k are the annihilation operators at A and B 
sublattices, respectively, 



H=0, 1,2,3 



(8) 



(To is the 2x2 unit matrix, <jj (j = 1, 2, and 3) are Pauli 
matrices and 

e (k) =~(e A + eg) - (t' c + ® cos( Vl • k) 

-(t' 6 +*fe)cos(v 2 -k) 

-(4 + ^)co S ((v 2 -v 1 )-k), (9) 

ei(k) =-t a cos(S a ■ k) - t b cos(S b ■ k) 

- t c cos(<j c ■ k) 

- t 3a cos(2<5 a ■ k) - t 3b cos(25 6 • k) 
-i 3c cos(2<5 c -k), (fO) 

e 2 (k) = - t a sm(5 a ■ k) - t b sin(d b ■ k) 

- t c sin(S c ■ k) 

- t 3a sin(-2<5 Q • k) - t 3b sin(-2<5 b • k) 
-t 3c sin(-2<5 c -k), (ff) 

6 3 (k) =l(e A - en) - (t' c - K) cos( Vl • k) 

- (t' b - €) cos(v 2 • k) 

-(^-Ocos((v 2 - Vl )-k). (12) 
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The energy is given by 



(a) 



£±(k) =e (k)± 



\ 



(13) 



j'=i 



The nearest-neighbor hoppings and the third-nearest- 
neighbor hoppings connect the sites in A and B sub- 
lattices, while the next-nearest hoppings connect sites 
within A sublattice or B sublattice. Therefore, the near- 
est hoppings and the third-nearest-neighbor hoppings ap- 
pear only in the off-diagonal elements (Eqs. (10) and 
(11)), and the next-nearest hoppings appear only in the 
diagonal elements (Eqs. (9) and (12)). 

There are no gaps at k = k*, when £i(k*) = e 2 (k*) = 
e 3 (k*) = 0. In general, each of the equations ej(k) = 
(i = 1, 2, and 3) gives the line in the plane of k x and k y . 
Therefore, it is only in rare case that three lines (cj(k) = 
(i = 1, 2, and 3)) intersect at the same point 21,22 . 
However, when the system has an inversion symmetry, 
i.e., A and B sublattices are equivalent, we get 



e 3 (k)=0, 



(14) 



and there exist Dirac points where two lines (ei(k) = 
and e 2 (k) = 0) intersect. Hereafter, we study the system 
with inversion symmetry, i.e., Eq. (14) is satisfied. 

As we have discussed, the next-nearest-neighbor hop- 
pings appear only in eo(k) and £3(k). Therefore, the 
next-nearest-neighbor hoppings do not affect the exis- 
tence and the location of the Dirac points in the system 
with inversion symmetry. They only affect the energy 
dispersion near the Dirac points, i.e. they may tilt a 
Dirac cone. Since we are interested in the location of the 
Dirac points and the topological phase transitions, we ig- 
nore the effects of the next-nearest-neighbor hoppings in 
this paper. 



III. TOPOLOGICAL NUMBER 

We define e(k) and <j>(k) by 

ei(k)=e(k) cos(0(k)), 
e 2 (k) = e(k) sim>(k)). 

Then the matrix given in Eq. (8) is written as 

S u = e(k) (cos(0(k))oi + sin(>(k))(72) . 

and it is diagonalized as 

U^SkUk = e(k)a 3 , 

where 

- 0(k) ■ -jr 

J7 k = e - l —< T3 e- l T 72 

/ . 4>(k) ■ £(k) \ 

1 / e 2 — e 2 \ 
= 71 { e i±S P j ' 



(15) 
(16) 



(17) 
(18) 



(19) 




(b) 




FIG. 2. (Color online) Lines for ei(k) = (thin blue lines) and 
£2(k) = (thick green lines). Parameters are t a = tb = t c = 1, 
tsa — tzb = tzc — in (a), and t a = fa = t c = 1, H a = 
tsb = t3c = 0.2 in (b). Orange circles show K, K' and their 
equivalent points. Small purple circles are Mi (i — 1,2 and 
3) and their equivalent points, and triangles are T and their 
equivalent points. Arrows show the vectors, (ei(k), 62 (k)). 
Dirac points appear at the intersection points of the blue lines 
and the green lines, which are K and K' in both figures. New 
green circular lines appear around k = 0, if ts a > t a /8 as 
shown in (b). The Dirac points are, however, only at the K 
and K', if ts a < t a /3. The topological numbers are +1 and 
— 1 for K and K', respectively. 



Therefore, the vectors, 



*±(k) 



(20) 



are the eigenvectors of the matrix with eigenvalues 
±e(k). 

In Fig. 2, we plot the lines given by £i(k) = (thin 
blue lines) and e 2 (k) = (thick green lines) for t a = tj, = 
t c = 1, ha = t 3b = t 3c = in (a) and t a = t b = t c = 1, 



ha = t 3b = t 3c — 0.2 in (b). The Dirac points are given 
by the intersection points of the blue lines and the green 
lines, which are at 

K = (0,y), (21) 

K' = (0,-f), (22) 

and their equivalent points related by the reciprocal lat- 
tice vectors. We also plot the vector (ei (k), e 2 (k)) in 
Fig. 2. The Dirac points are the cores of the vortex for 
the vector (ei(k), e 2 (k)). The topological number is +1 
and — 1 at K and K', respectively. Note that the config- 
urations of the arrows have a topological number +1, if 
they are clockwisely or counterclockwisely rotating flows 
around the core, or if all arrows are pointing out or into 
the cores, since </>(k) changes +27T when k goes counter- 
clockwisely around the core. The configurations of the 
arrows have a topological number —1, if they are point- 
ing out the core in ±x direction and pointing into the 
core in the ±y direction, for example. 

Since the eigenstates are given by </>(k) (Eq. (20)), the 
topological number ±1 is related to the Berry phase of 

±7T. 



IV. MERGING DIRAC POINTS AT M OR F IN 
THE BRILLOUIN ZONE 

When k* is a Dirac points, — k* is also a Dirac point if 
the system is invariant with respect to the time-reversal 
operation, which is always the case in this paper. Two 
Dirac points, which are exchanged by time-reversal oper- 
ation, can merge only at the points k* where k* and — k* 
are the same or are different by the reciprocal lattice vec- 
tors. These points are only four points in the Brillouin 
zone, Mi, M 2 , M 3 , and T, 



Mi 


V37T 

= ( 3 


(23) 


M 2 


= ( 2 f ,o), 


(24) 


M 3 




(25) 


r 


= (0,0). 


(26) 



From the equation 

e 1 (M 1 ) = e 2 (Mi)=0, (27) 

we obtain the condition for two Dirac points to merge at 
Mi as 

t a +t b -t c = t 3a + t 3b + t 3c : Mi . (28) 

Similarly we obtain the conditions for two Dirac points 
to merge at M 2 , M3, and T as 

-t a + h + t c = t 3a + t 3b + t 3c : M 2 , (29) 

ta-h + t c = t 3a +t 3b + t 3c :M 3 , (30) 
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(a) 



(b) 




FIG. 3. (Color online) Planes given by Eqs. (28) - (31) 
in the parameter space of t a , t b and t c in the case of (a) 
tsa = t 3b = t 3c = and (b) t 3a = t 3b = t 3c > 0. On the 
planes labeled by Mi, M2, M3 and V, two Dirac points merge 
at Mi, M2, M3, and T, respectively. The 3D space is divided 
into 14 regions, Ri - R4, R a , R b , Rc, Ri - Rt, R* a , R b , and R* 
when t 3a = t 3b — t 3c — 0. The regions, Rl, R2 etc. are given 
by the inversion of J?i, R2 etc.. When t 3a = t- ib = t 3c 7^ 0, 
a regular tetrahedral region Ro appears around the origin, as 
shown in the inset of (b). Topological phase transitions occur 
and the number of the Dirac points changes by two when 
parameters intersect the planes, Mi, M2, M 3 and T. 

and 

—t a — t b — t c = t 3a + t 3b + t 3c : r, (31) 

respectively. We plot the planes given by Eqs. (28) - (31) 
in the parameter space in t a , t b and t c in the cases of 
t 3a = t 3b — t 3c = and t 3a = t 3b — t 3c > in Fig. 3. 




(a) 



FIG. 4. (Color online) Phase diagram in the plane of tt/t a - 
t c jt a for the tight-binding electrons on honeycomb lattice 
with direction-dependent nearest-neighbor hoppings (t-^a = 
^36 = i'ic = 0). 



The phase diagram in tf,/t a — t c /t a plane (Fig. 4) can be 
obtained as the intersection of the phase diagram in 3D 
plot (Fig. 3) and the plane of constant t a , as shown in 
Fig. 5. 

If t 3a + t 3b + t 3c = 0, the t a — % — t c space is divided 
into 14 regions by four planes, Eqs. (28) - (31), as shown 
in Fig. 3 (a). There are two gapless Dirac points when 
parameters are in the regions labeled by R\, R2, R 3 and 
i?4 and their inversion regions, i?*, i?2, R 3 and R\ in 
Fig. 3. The regions, R\, i?2, R3 and i?4, correspond to 
the shaded regions in Fig. 4, where two-dimensional plot 
of the phase diagram is shown in tb/t a - t c /t a plane. On 
the planes and lines (Mi, M2, M 3 and T) in Fig. 3 and 
Fig. 4, topological phase transition occurs and two Dirac 
points merges at Mi, M2, M3 or T in the Brillouin zone. 
The phase diagram for t 3a + t 3 b + t 3c = is the same as 
obtained previously 7 . 

We show the 3D plot of the phase diagram in Fig. 3 (b) 
for t 3a + tsb + t 3c > 0. Although it is invisible in the 
main figure in Fig. 3 (b), there is a regular tetrahedron 
region Rq around the origin (0, 0, 0) as shown in the inset 
in Fig. 3 (b) . We plot the phase diagram in the plane of 
h/t a - tc/ta plane for t a = 1 and t 3a = t 3b = t 3c = 0.19 in 
Fig. 6. If the third-nearest-neighbor hoppings are finite, 
the number of the Dirac points can be larger than 2 and 
Dirac points can merge at the points other than Mi, M2, 
M3 and r. The topological phase transitions caused by 
these merging Dirac points are shown by the orange lines 
in Fig. 6, which we will discuss in Section VI. 

The boundary lines Mi and M 3 in tf,/t a — t c /t a plane 
move closer and the regions R\ and R2 becomes narrower, 
as t 3a +t 3 b+t 3 c becomes large, as seen in Fig. 5 (a). When 
*3o + ^36 + t 3c = t a , the boundaries (Mi and M 3 ) in the 
tb/ta ~ t c /t a plane overlap and the regions R\ and R 3 
disappear in the phase diagram in the tb/t a — t c /t a plane. 
If t 3a + t 3 b + t 3c > t a , Mi and M 3 lines are exchanged 
and the new region Rq appears, as seen in Fig. 5 (b). 



t 3 a=t 3 b=l3c=0-19 




(b) 



t 3a =t 3b =1 3c = 0- 4 




FIG. 5. (Color online) 3D plot of the phase diagrams for (a) 
tsa = t3b = t3 C = 0.19 and (b) 0.4 . The intersection with the 
plane t a = 1 is the phase diagram in Fig. 6 and Fig. 15 (b). 
Three Dirac points merge when parameters are at Pi, P2, or 
P3, and Ti, Ti, T2, and T3 are tricritical points, which are 
discussed in Section VI and Appendix. 



MOVING DIRAC POINTS WITH ONLY 
NEAREST NEIGHBOR HOPPINGS 



In this section, only nearest neighbor hoppings are 



considered (t' a = t' b = t' c = t" a - 
ha = t 3b = t 3c = 0). When \t a \ 
the "triangle inequality" , 



t'b = t" = and 
\tb\ and \t c \ satisfy 



N 
l*-l 



<M<M 
- t a - a + 



(32) 
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t = t, = L =0.19, t = 1 

3a 3b 3c a 




t a =l, t b =l,t c =(-2,0),(0,2) 




FIG. 6. (Color online) Phase diagrams in the ti/t a - t c /t a 
plane for t a = 1 and t^ a — t^b = <3c = 0.19. When the third- 
nearest-neighbor hoppings are finite, phase transition lines 
Mi, M2, M3 and Y are shifted and new transition lines appear, 
as discussed in Section VI and Appendix. The phase diagrams 
for other values of third-nearest-neighbor hoppings calculated 
numerically are show in Fig. 14 and Fig. 15. The numbers in 
the figures indicate the numbers of the Dirac points in the 
Brillouin zone. 



there are two Dirac points at k = ±k*. The position of 
the Dirac points is obtained by the equations 7,21 ' 22,26 



cosfvi • k* 



cos(v 2 • k*) = 



cos((vi - v 2 ) • k* 



t 2 


- t 2 - 


t 2 




2t a tb 




t 2 


- 1 2 - 


t 2 




2t c t a 




t 2 


- t 2 - 


t 2 

' c 


2tbt c 



(33) 
(34) 
(35) 



Note that the third equation (Eq. (35)) is satisfied if the 
first and the second equations (Eqs. (33) and (34)) arc 
satisfied and 



tb£ c sin(vi • k*) sin(v 2 • k*) > 0. 



(36) 



The condition, Eq. (32), is fulfilled in the regions R21 
R 3 , R 4 , R$, R%, R* 3 , and R\, in Fig. 3 (a). We plot 
the trajectories of the Dirac points as t c is changed for 
t a = tf, = 1 in Fig. 7. When t c = 1 two Dirac points 
are located at K and K' (up triangles and down trian- 
gles labeled by "1" in Fig. 7). The Dirac points move to 
Mi and the equivalent points as t c increases (blue broken 
lines and black lines in Fig. 7), and they merge at Mi 
and the equivalent points when t c = 2. Dirac points an- 
nihilate when t c > 2. Merging of two Dirac points causes 
a topological phase transition and this phase boundary 
is shown as Mi in Fig. 3 and Fig. 4. On the other hand, 



FIG. 7. (Color online) Locations and trajectories of the Dirac 
points as a function of t c with t a = tb = 1 being fixed. Only 
nearest neighbor hoppings are considered. Solid and Bro- 
ken lines are the trajectories of the Dirac points, which are 
converted by the space inversion each other. Up triangles 
and down triangles are the position of the Dirac points at 
t c indicated at the figure. When t c = 2 and t c = —2, two 
Dirac points merge at the Mi point (-^j: 71 ") and the F point 
(0,0), respectively. When t c = the system becomes one- 
dimensional. Therefore, the smallest up and down triangles 
are not the Dirac points, but Dirac points approach these 
points as t c — > 0. 



when t c is decreased, Dirac points move in the oppo- 
site direction. Dirac points approach to the midpoint of 
M 2 and M 3 , ±(M 2 + M 3 )/2 = ±(V3tt/2, -tt/2). When 
t c = 0, the system becomes one-dimensional and two 
bands intersect on the line. The smallest triangles with 
the label "0" in Fig. 7 should be considered as a singular 
points for the location of Dirac points as t c — > 0. As seen 
in Fig. 3 (a) and Fig. 4, t a = tb = 1 and t c = is the 
line (in Fig. 3 (a)) or the point (in Fig. 4) of intersection 
for the boundary planes or lines of M 2 and M3. When t c 
is changed to be negative, the parameters move into the 
region R4 in Fig. 3 (a) and in Fig. 4. Dirac points move 
on the cyan broken lines and the red lines in Fig. 7 and 
merge at T when t c = —2. 

Similarly, if we change t a with fixing tb — t c = 1 , Dirac 
points move on the vertical lines, i.e., the figure is rotated 
by 120 degree. In that case Dirac points merge at M 2 
and the equivalent points when t a = 2. 

In Fig. 8 we plot a set of trajectories as t c is changed 
(1*6 — t a \ < t c < tb + t a ) for several values of tb/t a > 0. 
Dirac points are located in the triangular regions around 
K and K', when the parameters t a , tb, and t c are in the 
R\ region in Fig. 3. Similarly, we obtain that the Dirac 
points move in other triangular regions in Fig. 9 when t ai 
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k 1% 



FIG. 8. (Color online) Locations of the Dirac points in 
the anisotropic honeycomb lattice with t a = 1, t& = 
0.1, 0.2, ■ • ■ 0.9, 1.0, 1.1, ■ • ■ 1.9, 2.0, 3.0, 4.0, • ■ • 9.0, 10.0. The 
lines are the trajectories of the Dirac points for \tb — t a \ < 
t c <t b +t a (region R± in Fig. 3) (a). 
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FIG. 10. (Color online) Lines for ei(k) = (thin blue lines) 
and £2(k) = (thick green lines) for t a = tb = to = 1, tz a = 
tsb — tic = 1/3. Green and Blue lines touch at Mi, M2 and 
Mi. 




FIG. 11. (Color online) Lines for ei(k) = (thin blue lines) 
and £2(k) = (thick green lines) for t a = tb = t c = 1, 
tsa ~ tsb = ts c = 0.4. Dirac points appear at K and K' and 
additional points (red circles, Di, D2 and D3). The topo- 
logical number is —1 and +1 for the additional three points 
around K and K', respectively. 



FIG. 9. (Color online) Divided regions of the moving Dirac 
points in the case of (t a ,tb, t c ) being in the regions of R\, R2, 
R3 and R4 in Fig. 3(a). 



VI. DIRAC POINTS IN THE SYSTEM WITH 
THIRD-NEAREST-NEIGHBOR HOPPINGS 



tb, and t c are in R 2 , R3 or R4 regions in Fig. 3. In Fig. 9 
we plot the regions for the Dirac points when parameters 
t a , tb and t c are in i?2, R3 and R4 in Fig. 3. 



When the third-nearest-neighbor hoppings are finite, 
more than two Dirac points are possible to exist. For 
simplicity we assume all hoppings, t a , tb, t c , t 3a , t 3 b, and 
t 3c are positive. Generalization to the negative hoppings 
is straightforward. 
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FIG. 12. (Color online) Lines for ei(k) = (thin blue lines) 
and £2(k) = (thick green lines) for t a = tb = t = 1, 
ts.a = t'ib = t-ic = 0.5. Four Dirac points around K and 
K' merge and the topological number becomes =f2 at K and 
K', respectively. 



v X \ H \ 

v \ \ \ \ 

2- \ \ \ 

K=D ft 





FIG. 13. (Color online) Lines for ei(k) = (thin blue lines) 
and £2(k) = (thick green lines) for t a = tb = tc = 1, t^ a = 
^36 = ts,c = 0.6. Dirac points locate at the intersection points 
of green lines and blue lines, i.e., K, K', Di, D2 and D3. 



A. Symmetric case (C§ v ) 



In this subsection we discuss the symmetric system, 



t — t a £ c , 

^3 = Ha = t3b = t^c 



(37) 
(38) 



Bena and Simon 15 have studied this model. They have 
shown that new Dirac points appear when £3 > t/3. As 
shown in Fig. 2, the blue lines, which is given by ei(k) = 
0, and the green lines, which is given by £2(k) = 0, have 
intersection points only at K and K', when £3 < t/3. 
Derivation is given in Appendix A. The blue lines and 
the green lines touch at Mj, M 2 and M 3 , when £ 3 = i/3 
(Fig. 10). We perform power series expansions of £i(k) 



and e 2 (k) about M2 and obtain 
d(k) ~+t(k x 



2V3tt 1 ,., 

k.. 



V3 , ( , 2V3tt 1, 2 



ea(k) ~ — t I k x - ; 




(39) 
(40) 



when t 3 — i/3. This means that merged Dirac point 
exists at M2 when £3 = t/3 and the energy is given by 



VMk)) 2 + Mk))2 




2%/3tt 



(41) 

in one direc- 



in other direction. 



i.e., the energy is proportional to \k x — k*\ 
tion and proportional to (k y — k*) 
The merged Dirac points, which have these energy dis- 
persion, are called semi-Dirac points 27 . When the semi- 
Dirac points exit, density of states is proportional to the 
square root of the energy 



D{ 



€) (X 



(42) 



and the energy of the Landau levels varies as 
B 2 / 3 [Ref. 28]. 

When £3 is larger than i/3, each touching point of 
blue and green lines changes to a pair of the intersec- 
tion points, i.e., each semi-Dirac point separate into two 
Dirac points with topological number ±1, and there exist 
eight Dirac points, as shown by Di, D 2 and D3 around 
K (= Do) in Fig. 11. As tj,jt becomes large, two Dirac 
points stay at K and K' and other Dirac points move 
from Mi, M 2 and M3 to K and K'. Three Dirac points 
moving to K for 1/3 < £ 3 /i < 1/2 (Di, D 2 and D 3 in 
Fig. 11) have the topological number —1 and three Dirac 
points moving to K' have topological number +1. When 
ts = t/2, four Dirac points merge at K and K'. In that 
case two blue lines and two green lines intersect at K and 
K', as shown in Fig. 12. The merged Dirac points at K 
and K' have topological number —2 and +2, respectively 
(1 x (+1) + 3 x (-1) = -2 at K, for example). The topo- 
logical number at K can be also obtained as follows. For 
(k x ,k y ) w K = (0, 47r/3), we write 



kcos(#), 



k y — — 7T = Ksin(#). 



When £3 = t/2, we obtain 



ei(k) fa -iK 2 cos(26>) 
8 

e 2 (k) ps --tn 2 sin(20). 
8 



(43) 
(44) 

(45) 
(46) 



Comparing these equations with Eqs. (15) and (16), we 
obtain 



0(k) = -29, 



(47) 
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and the topological number of the vector (ei(k), ei(k)) 
at K is —2. The energy at this point is obtained as 



(a) 



VMk))2 + ( £2 (k))2 



-t ki 




> (48) 



which has been obtained by Bena and Simon 15 . In this 
case the density of states is constant near e s» 0. 

When i 3 > f/2, eight Dirac points appear again. Two 
Dirac points are at K and K' having topological number 
+1 and —1, respectively. There exist three Dirac points 
with topological number —1 around K and three Dirac 
points with topological number +1 around K'. They 
form the equilateral triangles as in the case of t/3 < 
<3 < t/2, but the triangles are upside-down, as shown in 
Fig. 13 for t 3 = 0.6t. 



B. direction-dependent hoppings with 
third-nearest neighbor hoppings 

If the hoppings depend on the direction, the phase di- 
agram becomes much richer. Even if the hoppings be- 
tween the third-nearest-neighbor sites are small, there 
appear new phases as shown in Fig. 6. Two Dirac points 
can merge at the point other than Mi, M2, M3, and 
r. When a pair of the Dirac points merge at k*, which 
are neither Mi, M2, M3 nor T, in the Brillouin zone, 
other pair should merge at — k* simultaneously, since we 
study the system with time-reversal symmetry. There- 
fore, when the parameters are moved across these phase 
boundary lines (orange lines in Fig. 6 and curved lines 
in Fig. 14 and Fig. 15), the number of the Dirac points 
changes by four, i.e., 4-> 4, 2 -s-> 6 and 4 «->• 8. 

When the system has the reflection symmetry {C2V 
symmetry) we can obtain the analytical expressions for 
some special points (Pi, P2, P3, Qi, Q2, and Q3) in 
the phase diagram (see Appendix B). At Pi, P2 and P3 
in the phase diagram in tb/t a — t c /t a plane, three Dirac 
points merge. 

There exist the tricritical points, at which three phases 
with zero, two and four Dirac points terminate, in the 
parameter space in % and t c , which we label as Ti, T'i, 
T 2 and T 3 . 

In order to make the discussion simpler, we assume 
the third-nearest-neighbor hoppings are independent of 
the direction, i.e., we take £3 = t^ a — — t^ c . We plot 
Pi, P 2 , P 3 , Ti, T' 1; T 2 and T 3 as lines in the 3D plot in 
ta — % — t c space in Fig. 5. In the if, — t c plane at t a = 1, 
these points are obtained as 




and 



(1, 3*a 



(49) 
(50) 

(51) 
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FIG. 14. (Color online) Phase diagram in the tt-t c plane 
obtained numerically at t a = 1 and tz a = tzb = i3c = 0.25, 
0.3, and 1/3. The numbers of the Dirac points in the Brillouin 
zone is shown in the figure. 



Note that while Pi exists only for t 3 > 1/3, P2 and P 3 
exist for any value of £3 ^ 0. The points Qi, Q2 and Q3 
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FIG. 15. (Color online) Phase diagram in the tt-t c plane at 
t a = 1 and t 3a = t 3b = t 3c = 0.36, 0.4 and 0.4058. The 
numbers in the figure indicate the number of the Dirac points 
in the Brillouin zone. 



are given in Appendix B 2. We find that Qi exists when 



t 3 > 1/5, and Q2 and Q3 exist when * 3 > 1/4. 
The tricritical points are given as 



Ti 



l + 3* 3 , l + t 3 l-5t. 



l + 3t 3 



1 + 3t 3 l + *3 /l — 5i 



2 2 V 1 + 3< 3 

T' 1 = ((T 1 ) W) (T 1 ) ie ) J 

I 3* 3 , I -/; / l + 5t 3 
1 - 3t a ' 



-l + 3t 3 l-t 3 /l + 5t 



l-3t 3 



and 



((T 2 )„(T 2 ) X ) 



(52) 
(53) 

(54) 
(55) 



Note that Ti and T[ exist when < i 3 < 1/5, and T 2 
and T 3 exist when < * 3 < 1/3. 

The phase transition lines connections Pi, Qi and Ti 
are obtained numerically and shown in Figs. 14 and 15 for 
t 3 = 0.25, 0.3, 1/3, 0.36, 0.4 and 0.4058. The numbers of 
the Dirac points in the Brillouin zone changes by 2 when 
the phases are divided by the straight line and it changes 
by 4 when the phases are separated by the curved line. 

In Fig. 16 we show the regions where Pi and Qi locate. 
We plot the * 3 dependences of the coordinates of Pi, Qi 
and Ti in Fig. 17. 

At the tricritical points, four Dirac points with topo- 
logical numbers, +1, +1, —1 and — 1, merge at one of Mi 
in the Brillouin zone. 

The density of states are calculated in Appendix D. 
We obtain 



at Dirac point 

when two Dirac points merge 
(semi-Dirac point) 

when three Dirac points merge (56) 

Pi,P 2 ,P 3 

at tricritical point 

Ti.Ti.Ta.Ta. 



Note that although the topological number of the merged 
point of three Dirac points is the same as that of the 
ordinary Dirac point, energy dependences of the density 
of states are different. 

In Fig. 18 we show how Dirac points move when t a is 
changed when other parameters are taken as = t c = 1 
and * 3a = < 3 b = < 3c = 0.4. The 3D plots of the energy are 
given in Fig. 19 and Fig. 20. When t a = 1, Dirac points 
are locates at K and apexes of the regular triangle (black 
circles in Fig. 18. See also Fig. 11.). As t a becomes 
larger, two Dirac points on the k y axis comes closer and 
merge when t a = 41/40 = 1.025 (this set of parameters 
corresponds to Q2 or Q 3 in Fig. 15 (b) with exchanging t a 
and t\y (or * c )). These Dirac points disappear when t a > 
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FIG. 16. (Color online) Locations of the parameters for three 
Dirac points to merge (Pi, P2 and P3) and for two Dirac 
points to merge in the axisymmetric case (Qi, Q2 and Q3). 
When £3 = 0.5 all points are at (1,1). They do not exist in 
the shaded region (0 < t 3 < 1/5 for Qi, < t 3 < 1/4 for Q 2 
and Q3, and < t 3 < 1/3 for Pi). The regions labeled as 
R a , Rb, Rc, Ro, R2, Ri and R 3 , are those shown in Fig. 3. 
The vertical lines labeled as Mi, M 2 and M3 show that Pi 
and Qi (i — 1, 2, and 3) are on the boundaries (Mi, M2 and 
M3, respectively). 
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FIG. 17. (Color online) Coordinates of the special points in 
the phase diagram in tb/t a - t c /t a plane as a function of £3 
(£3 = t 3a — = t 3c and t a = 1). At Pi, P2 and P3 
(Eqs. (49), (50), and (51)), three Dirac points merge, and at 
Qi, Q 2 and Q 3 (Eqs. (B16), (B17), and (B18)), two Dirac 
points merge on the symmetric line (k x — 0) in the Brillouin 
zone. 



1.025. Other Dirac points move to Mi and M3 and they 
merge when t a = 1.2. When t a becomes smaller, three of 
the four Dirac points near K merge when t a = 0.944 (this 
set of parameters corresponds to P2 or P3 in Fig. 15 (b), 
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FIG. 18. (Color online) The moving Dirac points by changing 
t a with fixed t b = t c = 1 and t 3a = t- ib = tic = 0.4. 



see also Fig. 21) and one Dirac points remains after three 
Dirac points merge (t a < 0.944). Another one of the four 
Dirac points move to M2 and merge when t a = 0.8. 



VII. CONCLUSIONS 

We study the Dirac points of the electrons on the hon- 
eycomb lattice with up-to third-nearest neighbor hop- 
pings. Since the next-nearest-neighbor hoppings do not 
affect the location of the Dirac points, we study the 
location of the Dirac points and the phase diagram in 
the parameter space in the direction-dependent nearest- 
neighbor hoppings (t Q , tb and t c ) and the third-nearest- 
neighbor hoppings (t^a, and tj, c ). Dirac points are 
obtained as the intersection points of two kinds of lines, 
which are given by ei(k) = and £2(k) = 0. Dirac 
points move from K and K', where Dirac points exist 
if hoppings are independent of the direction. We ob- 
tain the trajectories of the Dirac points as the direction- 
dependent hoppings are changed. The topological phase 
transitions occur when Dirac points merge and disap- 
pear. We obtain the phase diagram in the parameter 
space of the direction-dependent hoppings. In each phase 
there are 0, 2, 4, 6, or 8 Dirac points, half of which 
have the topological number +1 and the other half have 
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FIG. 19. (Color online) 3D plot of the energy as a function 
of the wave number for the parameters in Fig. 18 (t a > 1). 




FIG. 20. (Color online) 3D plot of the energy as a function 
of the wave number for the parameters in Fig. 18 (t a < 1). 



the topological number —1. The number of the Dirac 
points in the Brillouin zone changes by 2 when two Dirac 
points merge at Mi, M 2 , M 3 , or T, and it changes by 
4 when two pairs of Dirac points merge at other points 
in the Brillouin zone, ±k*. When parameters are at Pi, 
P2 or P3, three Dirac points with the topological num- 
ber +1, —1, and —1 (or —1, +1 and +1) merge. Four 
Dirac points with topological number +1, —1, —1, and 
— 1 (—1, +1, +1, and +1) merge at K (K') only when 
t a = tf, = t c = 2t^ a = 2i3b = 2<3 C . We also obtain that 
there are tricritical points Ti, T'i, T2 and T3, where 
three phases with 0, 2 and 4 Dirac points terminate. 
There are no parameter regions where two Dirac points 
with the same topological number merge. 

The density of states D(e) for e « is proportional 
to |e|, y/e, e 1 / 3 , e 1 / 4 , and constant at the Dirac point, 
the merged point of two Dirac points (with topological 
number ±1 and ^1), the merged point of three Dirac 
points (two of them have the topological number ±1 and 
the other has the topological number =Fl), the tricritical 
point, and the merged point of four Dirac points (three 



of then have the topological number ±1 and the other 
has the topological number =f1), respectively. 

In order to realize these topological phase transitions, 
only small third-nearest-neighbor hoppings are neces- 
sary, if the anisotropy is sufficiently large. The third- 
nearest-neighbor hoppings in monolayer graphene are es- 
timated to be finite (about 0.1 t) 20 . On the other hand, 
if there are large third-nearest-neighbor hoppings, small 
anisotropy is enough to realize these topological phase 
transitions. The strained bilayer graphene is considered 
to have large effective third-nearest-neighbor hoppings 14 . 

Therefore, these topological phase transitions can 
be realized in the ultracold atoms on the optical lat- 
tice, strained monolayer graphene or strained bilayer 
graphene. 



13 



Appendix A: Dirac points in a symmetric system 

We study the system with Cq v symmetry in this Ap- 
pendix, i.e., we set 



t — t a t b t c 



(Al) 



and 



h = ha = hb = he- (A2) 
Then ei(k) and e 2 (k) (Eqs. (10) and (11)) are written as 



ci(k) = -t 



h 



V^3 1 
cos( — k x ) + 2 cos( — k x ) cos(-k y ) 

3 z 

, 2\/3 . . \/3 . , . 

cos(— — k x ) + 2 cos( — k x ) cos{ky) 



(A3) 



and 



e 2 (k) 



is 



. . \/3 . . , \/3 . . 1 . 
sm( — k x ) - 2sm( — k x )cos(-k y ) 

. ,2\f3 . . ,\/3, % /, x 
sin(— — k x ) - 2 sm(—k x ) cos(k y ) 



(A4) 



In Figures 2, 10, 11, and 12, we plot the lines for ei(k) = 
(thin blue lines) and for e2(k) = (thick green lines). 
The massless Dirac points are given by the intersection 
points of the blue lines and the green lines. We obtain 
62 (k) = if k x — from Eq. (A4). From the symmetry 
of the system, we obtain that £2(k) = on the thick 
green lines in Fig. 2 (a). When £3 = there are no other 
solutions of e 2 (k) = 0. 

In order to find other solutions of e 2 (k) = in the case 
of t 3 7^ 0, we expand e 2 (k) around k x = and obtain 

e 2 (k) w ^ ^-t + t cos ^ + 2t 3 - 2t 3 cosk^j k x 

+ 0{kl). (A5) 

The intersection points of green lines (e 2 (k) = 0) and k y 
axis are obtained by the equation 



fc (2) 

-t + t cos -|- + 2t 3 - 2t 3 cos k y 2) 



We obtain 



,(2) 



cos ■ 



t± \t-8t 3 \ 
8^ ' 



0. (A6) 



(A7) 



(2) 

This equation has the solution k y — 2mr with integers 
n, and other solutions exist if t 3 > t/8. As seen in Fig. 2, 



10, 11, and 12, there are circular-like green line around 
k = 0. Note that all T points are not equivalent if we 
consider ei (k) and e 2 (k) separately, while they are equiv- 
alent if we consider -\/(ei(k)) 2 + (e 2 (k)) 2 . 

As seen in Fig. 2(b) and Fig. 10, if t 3 < t/3, Dirac 
points are only at K and K', since green lines and blue 
lines intersect only at K and K'. In order to examine 
the number of Dirac points, we also expand ei(k) with 
respect to k x around k x = 0. Then we obtain 



(A8) 



ei(k) rj -t ^l + 2cos^-J -t 3 (l + 2cosfc y ) 
+ 0(k x ). 

We obtain from €i(k) = at k x = that 



cos ■ 



-t±\t- 2t 3 
4t 3 



(A9) 



which gives the intersection points of blue lines (ei(k) = 
0) and k y axis. We obtain that the solution of Eq. (A9) is 
only ky = ±47r/3 + 4n7r (K and K', and their equivalent 
points) if i 3 < t/3. If t 3 > t/3, there exist other solutions 
of Eq. (A9), which arc the intersection points of blue 
ellipse and k y axis, which we call Di as shown in Fig. 11 
for t 3 = 0.4t. Because of the symmetry the new Dirac 
points for t 3 > t/3 make an equilateral triangles around 
K and K', as shown in Fig. 11 (Di, D 2 and D 3 ). If 
i 3 = t/2, four Dirac points merge at K and K' as shown 
in Fig. 12. 



Appendix B: Dirac points in an axisymmetric 
system (C2v) 

In this appendix we study the system, which has C 2t , 
symmetry, i.e. 



tb t c — tbc, 
t 3 b = he = hbc- 



(Bl) 
(B2) 



In this case the system has n rotational symmetry and 
the reflectional symmetry with respect to the x-axis and 
the vertical line x — \/3/6. We obtain, 

v^3 v^3 1 

d(k) = - t a cos( — fc x ) - 2t bc cos( — k x ) cos(-k y ) 

. 2\f3 . . \p3 . . . 

- ha cos(— — k x ) - 2t 3fcc c °s( — k x ) cos(ky), 

(B3) 

and 

e 2 (k) =-t a sin(— fc x ) + 2t bc sin(—k x ) cos(-k y ) 

2\[3 . . . \[3 . . . 

+ ha sm(—^-k x ) - 2t 3bc sm( — k x ) cos(k y ). 



(B4) 
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It is obtained that 62 (k) = when k x = and some of 
the Dirac points are on the k y axis, as in the Cq v case 
(Appendix A). As in Appendix A, the intersection points 
of the green lines (e2(k) = 0) and the k y axis are obtained 
by 



x/3 / k {2) \ 

— \-t a + t hc cos -V— + 2t 3a - 2t 3bc cos k y 2) = 0. 



We obtain 



(B5) 



4 2) _ t bc ± y/^ - 16t a t 3bc + 32t 3a t 3bc + 32t 



cos 



Mir 



8f 



Mir 



(B6) 

The intersection points of the blue lines (ei(k) = 0) and 
the k y axis are obtained by 



k 



(i) 



t a - 2t bc cos - t 3a ~ 2t 3bc cos fcW = 0. (B7) 



(a) 



(b) 



_ 

K 



V V V V V 

\\\\\ 




K 



From this equation we obtain 



k y 1} _ -t bc ± y/t 2 bc - 4t a t 3bc ~ 4t 3a t 3bc + 8t 



cos 



Mk- 



4t 3bc 



1. three Dirac points merge 



(B8) 



If the solutions of Eq. (B6) and Eq. (B8) exist and 
they are the same, three Dirac points (D , D 2 and D 3 
in Fig. 11) merge at the k y axis as shown in Fig. 21. We 
obtain that three Dirac points merge on the k y axis when 
the equation 



(ta — ha — 2t 3bc )t bc + 4tl a t 3bc = 



(B9) 



is satisfied. Thus, we obtain that three Dirac points 
merge on the k y axis if (t a ,t b (= t bc ),t c (— t bc )) is given 
by 



(j^a •> 1 ^c) ( 1 



V~ta + t%a + 2t3bc V~^a + ^3a + ^3bc 

(BIO) 




Drf=D 2 -I^3 ^ ^ 




FIG. 21. (Color online) Lines foe ei(k) = (thin blue lines) 
and £2(k) = (thick green lines) for (a) t a = 1, t b = t a = 
1.13137, t Sa = t 3b = t Sc = 0.4, and (b) t a = 0.944, t b = t c = 1, 
t-ia ~ t3b ~ tzc = 0.4. In both cases three Dirac points merge 
on the k y axis near K and K'. The set of parameters in (a) is 
Pi in Fig. 15 (b), and the set of parameters in (b) corresponds 
to P2 and P3 in Fig. 15 (b) with exchanging t a and t b . With 
the anisotropy of t a and t b — t c , three Dirac points (Do, D2 
and D3) move from those in Fig. 11 and merge. 



(ta, t b , t c ) — ( t 3a + 2t 3bc — 4 3 ° 2 3bC ,tfe e , t bc ) . (Bll) 



"be 

Eq. (B10) is real when t 3a + 2t 3bc > t a . 

In Fig. 6, Fig. 14 and Fig. 15, we plot the phase 
diagram in the t b - t c plane by taking t a = 1 and 
t 3a = t 3bc = t 3 . We find three Dirac points merge when 
tb and t c are at the Pi, P2, or P3, which are given in 
Eqs. (49), (50) and (51). In order to derive Pi we take 
t a = 1 in Eq. (B10), and in order to derive P2 (P3) we 
take t bc = 1 in Eq. (Bll) and exchange t a and t c (t a and 
t b ). For example, if we take t 3a = t 3bc = 4/10, we obtain 
three Dirac points merge on the k y axis if (t a ,t b ,t c ) — 



(1,4^2/5,4^2/5) « (1,1.13137,1.13137), which corre- 
sponds to Pi in Fig. 15 (b), and (t a ,t b ,t c ) = (0.944, 1, 1), 
which corresponds to P2 and P3 in Fig. 15 (b). 

In Fig. 21 we show these two cases. At these param- 
eters three Dirac points, Do, D2 and D3, merge on the 
k y axis, as shown in Fig. 21. When parameters are those 
of Fig. 21 (a), which correspond to Pi in Fig. 15 (b), 
Di does not exist, while Di exists when parameters are 
those of Fig. 21 (b), which correspond to P2 and P3 in 
Fig. 15 (b). This difference can be understood by noting 
the following fact in Fig. 15 (b): When parameters move 
from the point t a = t b = t c = 1 to Pi, they intersect the 
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t„=t, h =t,»=0.3 



> \ X \ X > 

D<j=Dk \ \ \ 




FIG. 22. (Color online) Lines foe ei(k) = (thin blue lines) 
and t2 (k) = (thick green lines) for t a = 1, tb = t c = 0.980 
and t3a = t-ib = tsc = 0.4. This set of parameters is Qi in 
Fig. 6 (c). 



line M2, which shows that Dirac points merge and anni- 
hilate at M2 in the Brillouin zone when they cross the 
line M<2- On the other hand parameters do not intersect 
any lines labeled M\, M2 and M3 when they move from 
the point t a = t b = t c = 1 to P 2 or P 3 . 

When parameters are changed through these critical 
values, three Dirac points merges and one Dirac point 
survives. The topological number (Berry phase) of the 
Dirac point at Do is +1, while these at Di, D2 and D3 
are —1. As a result, the merged point of three Dirac 
points has the topological number —1. However, the 
merged point is not a simple Dirac point. The energy 
changes linear in one direction in the wave number and 
cubic in other direction, and the density of states is pro- 
portional to e 1 / 3 , as we show in Appendix D. 



2. two Dirac points merge 

Two Dirac points merge on the k y axis when two solu- 
tions of Eq. (B8) become the same, which happens when 

tbc ~ ^Aatzbc — ^tiat^bc + 8t^ bc = 0, 



and 



\tbc\ < 4|t 3 6c 



(B12) 



(B13) 



Thus, we obtain that two Dirac points merge on the k y 
axis if (t a ,tb(= tbc),t c [= t bc )) is given by 

t 2 

(t a ,tb, t c ) = (— < 3a + 2^3b c , tbc, he), (B14) 

4t3&c 




t =1 t =1.45 t b =1.35 




0.10 0.05 0.05 0.10 



0.10 0.05 0.05 0.10 t 



FIG. 23. (Color online) Example for a pair of Dirac points 
merge at M. Lines of ei(k) = (thin blue line) and £2(k) = 
(thick green line). 



When t^a = t^ c = and t a = 1, we obtain the special 
points in the phase diagram in t b — t c plane as, 

Qi = (203(1 -* 3 ), 203(1 -t 3 )) :(t 3 > 1/5), 

(B16) 

:(*3 > 1/4). 

(BIT) 



Q 2 =(t 3 + ^,i 



and 



Q3=(l,*3 + i - 



:(*3 > 1/4). 



(B18) 



For example, if we take < 3a = t 3 h c = 4/10, 
we obtain two Dirac points merge at (t a ,t{,,t c ) — 
(l,2V6/5,2V6/5) « (1,0.980,0.980), which corresponds 
to Qi in Fig. 15 (b), and at {t a ,t b ,t c ) = (41/40,1,1) = 
(1.025,1,1), which corresponds to Q 2 and Q 3 in 
Fig. 15 (b). In these values of parameters, two Dirac 
points Do and Di merge on the k y axis as shown in 
Fig. 22. 

We obtain the phase diagram numerically in Figures 
14 and 15. For example, in Fig. 15 (b), the triangular- 
like region enclosed by the line, P1-S1-Q3-P2-Q1-P3-Q2- 
S2-P1, is the parameter region of tb and t c where the 
additional Dirac points exist when t a = 1, and t 3a = 
t?,b = t-3c = 0.4. Therefore, the apexes of the triangular- 
like region in Fig. 15 (b) are Pi = (1.131, 1.131), P 2 = 
(0.994,1) and P 3 = (1,0.994), and the line passes on 
Qi = (0.980, 0.980), Q 2 = (1.025, 1) and Q 3 = (1, 1.025). 



Appendix C: tricritical points 



(fa, h, tc] 

- {ta, 2\/ (t a + t3a ~ 2t3bc)t3bc, 2\/ (t a + t 3a — 2t3bc)t3bcj ■ 

(B15) 



In this appendix we give the analytical expression for 
the tricritical points. When i 3a + 2t 3 b c < t a , Pi does 
not exist in the phase phase diagram in the tb/t a - t c /t a 
plane. In this case the lines P2-Q3 and P3-Q2, on which 



16 




o.io 0.05 o 0.05 o.io * 0.10 0.05 0.05 0.10 71 



expand £i(k) and £2(k) around k = (0, 2ir) and obtain 

ei(k) ps a(fc y - 2vr) 2 + /3(fc y - 27r)fc a + 7 *£, (C2) 
e 2 (k) «<5i(fc y -2 7 r)- ( 5 2 fc ;c , (C3) 

where 



a = 


1 


+ 3t 3b + M 3c ) , 


(CM) 


= 




+ t 3a — 3^3b + 5t 3c ) , 


(C5) 


7 = 


1 


" t 3 b + ^3c) , 


(C6) 


£1 = 


- (t„ - 2t c + 


t 3 a + 3t3b — i3 c ) , 


(C7) 



and 



FIG. 24. (Color online) Example for a pair of Dirac points 
merge at M. Lines of ei(k) = (thin blue line) and £2(k) = 
(thick green line). 
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FIG. 25. (Color online) Example for a pair of Dirac points 
merge at M. Lines of ei(k) = (thin blue line) and £2(k) = 
(thick green line). 



two Dirac points merge, terminate tangentially on the M\ 
and M3 line at the tricritical points, which we call T3 and 
T 2 , respectively (see Fig. 14). When 6i3f, c — t 3a <t a , Qi 
does not exist, and tricritical points (Ti and T^) appear 
on the M 2 line. In the following we give the analytical 
expressions of Tj., T^, T2 and T3. 

We study the case 



tb — t a — t c + t 3a + t 3 b + tic 



(CI) 



8 2 = ^-{ta-ha+tsb+tsc)- (C8) 

The line given by e 2 (k) = goes through (0, 2tt) and it 
is obtained as 



hy 2tt ~ k x 
Si 



(C9) 



On the other hand there are no lines near (0, 27r) given 
by ei(k) = 0, if (3 2 — < 0. We show the example in 
Fig. 23 (t a = 1.45, t b = 1.35, t c = 1, and t 3a = t 3b = 
t 3c = 0.3). If tb is changed smaller, the semi-Dirac point 
becomes a a pair of Dirac points near k = (0, 2tt) as 
shown in the left figure in Fig. 23. 

If j3 2 — 4a7 > 0, ei(k) = gives two lines, 



fc„ — 2-7T 



-p± V/? 2 -4o 7 . 
2a 



(C10) 



as shown in Figs. 24 and Figs. 25 ((t a ,t b , t c ) = (1.6, 1.5, 1) 
and (1.8, 1.7, 1), respectively, and t 3a = tab = He = 0.3). 

As seen in the left figures in Fig. 24, the merged Dirac 
points at k = (0, 2-k) disappear or become two Dirac 
points, when we change tb larger or smaller, respectively. 
Since there are other Dirac points in Fig. 24, this topolog- 
ical phase transition separate the phases with two Dirac 
points and four Dirac points. On the other hand, as seen 
in the left figures in Fig. 25, the merged Dirac points at 
k = (0, 2ir) become two Dirac points or disappear, when 
we change tb larger or smaller, respectively. Since there 
are no other Dirac points in Fig. 25, the topological phase 
transition is the transition between two phases with two 
Dirac points and zero Dirac points in this parameters. 
The tricritical point, on which three phases with zero, 
two and four Dirac points terminate in the parameter 
space, happens when one of the lines ei(k) = coincides 
with the line e 2 (k) = 0, 



Then merged Dirac point (semi-Dirac point) exists at M 2 
and the equivalent points, for example, k = (0,27r). We 
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Tricritical points Ti and can be obtained by solving 
Eq. (Cll) with respect to t c as a function of t a , t 3a , hb 
and he- We obtain 

7 t a + t 3a + t 3b + t 3c 2(hb — hc)ha 

tc± = ^ 1- 



± 



2 t a + t 3a + t 3b + t 3c 

ta - ha + t 3 b + h c 



x ta + ha — 3(t 3 b + h c ) f 2(t 3 b — h c ) 



t a + t 3a + hb + t 3c \t a + ha + hb + h c 

(C12) 

Since Ti and are on the line M 2 in the phase diagram 
in tb — t c plane with t a — 1, we obtain from Eq. (C12) 
and taking t a = 1 

(Ti) x = (1 — t c - + t 3 a + ^36 + t 3c )\ ta= i 

_ 1 + t 3a + hb + he 2(t 3o — t 3c )t 3 a 



+ 



2 1 + ha + hb + hc 

1 — ha + hb + t 3c 



„ / l+*3o-3(t 3 fc + *3c) i ( 2(t 3 b-he) 



1 + t 3a + hb + t 3c \ 1 + i3 a + £36 + t 3c J ' 

(C13) 



(TO, - (i c -)k=i 



1 + ha + ^3fc + he 2(t 3 ft — i3c)^3a 



2 1 + h a + t 3 b + he 

1 - ha + hb + t 3c 



x J l + ha -3(hb + he) | f 2(t 3fc -t 3c ) 



V 1 + t3a + *36 + t 3c \ 1 + ^3o + ^36 + ^3c 

(C14) 

(T'i)x = (1 - t c+ + t 3 a + hb + hc)\t a =l 

_ 1 + ^3a + % + he _ 2(t 3 fc - t 3c )t 3 g 

2 1 + *3a + hb + *3c 

1 — ha + ^36 + *3c 



x J l + ha -3(t 3 b+he) | f 2(t 3fc -t 3c ) 



1 + t 3a + hb + he \ 1 + t 3a + t 3 b + h 

(C15) 



We obtain Eqs. (52) and (53) by taking t 3a = hb = h c — 
t 3 in Eqs. (C13), (C14), (C15) and (C16). 

The tricritical points T2 and T3 are obtained as fol- 
lows. By solving Eq. (Cll) with respect to t a , we obtain 
the critical value of t a as a function of t c , ha, hb and t 3c 
as 

t c — t 3a — hb — t 3c 2(t 3a — hb)hc 



t a = 
+ 



2 t c — t 3a — hb — t 3c 

t c - ha — hb + hc 



ltc + 3t 3a + 3t 3b -t 



3c , / %(ha — hb) 

A \l h 

y t c — t 3a — hb — he \t c — ha — hb — h 

(C17) 

The tricritical point T2 is on the line M 3 in the phase 
diagram in h — t c plane, i.e., the merged Dirac points are 
located at M 3 instead of M 2 . Since the above calculation 
has been done for M 2 , we change t a — > h — > t c — > t a and 
ha — > hb —¥ t 3c —¥ t 3a , cyclically, in order to obtain T 2 . 
Then taking t a — 1, we obtain 

,m \ 1 — h a — t 3 b — he 2(t 3 b — t 3c )t 3a 

( i -2)x = ^ h 



+ 



2 1 — ha — hb — h c 

1 + t 3a - hb - he 



x ^ j 1 - t 3a + 3t 3o + 'it 3 c ( 2(hb — hc) 



1 — ha — hb — hc \ 1 — ^3o _ % — *3 

(C18) 



and 



(T 2 )j/ — (T 2 );e — 1 + t 3 a + ^3o + ^3c 

1 ^ t 3a — hb — t 3c 2(hb — hc)ha 



+ 



2 1 — t 3a — hb — hc 

1 + t 3a — hb — hc 



1 1 - ha + 3t 3 fc + 3t 3 



2(*3b ~ hc) 



V 1 ~ ha — hb —t 3c \ 1 — t 3a — t 3 b — h 

(C19) 

We obtain Eq. (54) by taking t 3 a — hb = h c = t 3 m 
Eqs. (C18) and (C19). 

We can obtain T 3 by changing t 3a -> h b -> t 3c -> t 3o , 
cyclically, in Eqs. (54) and (54). Then we obtain Eq. (55) 
similarly. 



and 

(Ti), - (? c+ )k=i 

*o + *3a + *3b + *3c 2(*3fo — t3c)*3a 



+ 



2 *a + h a + hb + he 

ta — t 3a + t 3 b + t 3c 



' t a + hg - 3(t 3 b + he) ( 2(hb - hc) 

ta + h a + t 3b + hc \t a + ha + t 3b + t 3c 

(C16) 



Appendix D: Density of states 

We study the density of states D(e) for |e| -c 1 in this 
Appendix. 



1. Density of states due to Dirac points 

First we study the Density of states due to Dirac point. 
A Dirac point (k*) is given by the intersection points 
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of two lines, ei(k) = and 62 (k) = 0, in the Brillouin 
zone. Two lines can be approximated as straight lines 
near the Dirac points. We take power series expansions 
of ei(k) and £2(k) about k*. By rotating the k x - k y axes 
appropriately and taking k* = for simplicity we can 
write 



ei(k) ~ -^ C ( k * + uk y), 
£2(k) w -^=C'(k x - uk y ), 



(Dl) 
(D2) 



where C(> 0), C'(> 0) and u(> 0) are constants. The 
density of states due to this Dirac point is calculated as 



D(e) 



1BZ if 



dkxdkyS 



( ei (k)) 2 + (£ 2 (k)) 2 



1 



S BZ CC'u' 



(D3) 



where Sbz is the area of the Brillouin zone. We obtain 
that the density of states is proportional to |e|. 



2. Density of states due to merged Dirac points 

The merged Dirac points at k = k* are classified into 
four types; 

(a) Two lines which are given by t\ (k) = and ti (k) = 
touch in the order n at k = k* . 

(b) One of the lines makes a loop and it shrinks into a 
point at k = k*. 

(c) The equation £2(k) = gives two intersecting lines 
at k = k* and the line given by the equation 
ei(k) = touches with one of the line in the or- 
der mat k = k* . 



(d) Both equations, ei(k) = and e 2 (k) 
lines intersecting at k = k* . 



give two 



There may be other possibilities, for example more than 
two lines are given by the equation £2(k) = in type (c) 
and (d), or two lines given by £2(k) = touch each other 
with order £(> 2) at k = k* in type (c). We neglect 
these possibilities for simplicity, since these cases are not 
realized in our present study. 

When t a = % = t c = 1 and ts a = = t$ c = 1/3, two 
Dirac points merge at Mi, M2 and M3 (see Fig. 10). At 
k = (2\/37r/3, 0), two lines touch quadratically (n = 2), 
which is type (a). At k = (0,7r), which is the equivalent 
point to k = (2a/37t/3, 0), the line given by the equation 
ei(k) = shrinks into a point (compare with Fig. 11), 
which is type (b). Although the merged Dirac point at 
k = (0,7r) looks different from that at k = (2-^/3^/3,0), 
the energy around these points ^/(e^k)) 2 + (e 2 (k)) 2 is 
the same. Type (c) of the merged Dirac points can be 
seen in Fig. 22, where two Dirac points, D and Di, 



merge on the k y axis when t a = 1, tb = t c = 0.980, and 
ha = hb = He = 0.4. 

When the parameters are at the tricritical points in 
the phase diagram (Ti, T^, T 2 and T 3 ), ei(k) = gives 
two lines and one of the lines touches cubically with the 
line given by 62 (k) = (see the point k = (0,27r) in 
Fig. 24 and Fig. 25). This is type (c) with cubic touching 
(m = 3). Four Dirac points with the topological numbers 
+1, +1, —1 and —1 merge at this point. In the same 
parameters for the tricritical point, type (a) with n = 4 
is realized at k = (2-^/3^/3,0), (not shown). 

First we study the density states for the merged Dirac 
points of type (a). If two lines touch in the order n, we 
perform the affine transformation, i.e. rotation, transla- 
tion, and stretching in the momentum space, and we can 
write 



e{ a) (k)^^=(k x +u n k^), 
e 2 •'(k) « —j={k x — u n ky), 



(D4) 
(D5) 



where u n (> 0) is a constant. The topological number at 
this point is ±1 if n is odd and if n is even. In this case 
we obtain 



^(eW(k)) 2 +(e^(k)) 2 = v ^ 
The density of states is obtained as 



1 20Fr( 



2n+l ' 

2n > 



S'bz K)^r(^) 



(D6) 



(D7) 



where S' BZ is the area of the Brillouin zone after the 
affine transformation. We obtain the Dirac point by tak- 
ing n = 1. When two Dirac points merge, we obtain a 
semi-Dirac point (energy depend linearly on the one di- 
rection of the momentum and quadratically in the other 
direction), which is obtained by taking n = 2. If three 
Dirac points merge, n — 3 is realized, as seen in the 
merged points near k = (2-^/3^/3,277/3) in Fig. 21. At 
the tricritical point in the phase diagram in the param- 
eter space, n = 4 is realized (at k ~ (2-^71-/3,0) for the 
parameters in Fig. 24 and Fig. 25, although that point is 
not shown in these figures. 

Next we study the density of states due to the merged 
Dirac point of type (b). Type (b) of the merged Dirac 
points can be described after the affine transformation as 



ei h \k)^C"(kl + kl), 



(D8) 
(D9) 



where C"(> 0) is a constant. In this simplification 



e^(k) = gives a straight line and = gives the 

point k = (0,0). We obtain ^(e^k)) 2 + (4 &) (k)) 2 = 
\k x \ +0(fc 2 ) in the k x direction and C'ky in the k y direc- 
tion. Therefore, Eqs. (D8) and (D9) describe semi-Dirac 



(*>), 
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points and the density of states is proportional to yfe. 
Indeed, the density of states is calculated as 



£>< 6 >(e)oc f 

J Inn 



2|e| 



2JU f 1 



zdk x 



-Ax 



(D10) 



where integral should be performed where the integrand 
is real and we have neglected higher order terms in e. 



The merged Dirac point of type (c) is described by 



e\ ''(k) ~ k x , 

e 2 ^(k) ~ ky(k x — u k' 



(Dll) 
(D12) 



In this case (k) = gives a straight line and e^' (k) = 
gives an intersecting straight line and a curve touching 
to the first line at k = (0, 0) in the mth order (if to 

I . e^(k) = gives two intersecting lines). We obtain 

( ei c) (k)) 2 + (4 c) (k)) 2 = \k x \ in the k x direction and 
u'\k y \ m+1 in the k y direction. Therefore, type (c) with 
to is similar to type (a) with n = m + 1. We obtain the 
density of states due to this merged points as 



D(e) 



(c) 



cx e "•+! , 



(D13) 



The simplest case of the type (d) is realized in the case 
t a = % = t c = 1 and t 3a = t 3 h = t 3c = 1/2 as discussed 
in Section VIA. In that case the topological number is 
two and the density of states is independent of energy for 
|e|«0. 
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